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ABSTRACT 

Investigating the molecular gas in the inner regions of protoplanetary disks provides insight into how 
the molecular disk environment changes during the transition from primordial to debris disk systems. 
We conduct a small survey of molecular hydrogen (H 2 ) fluorescent emission, using 14 well-studied 
Classical T Tauri stars at two distinct dust disk evolutionary stages, to explore how the structure of 
the inner molecular disk changes as the optically thick warm dust dissipates. We simulate the observed 
HI-Lyman a-pumped H 2 disk fluorescence by creating a 2D radiative transfer model that describes 
the radial distributions of H 2 emission in the disk atmosphere and compare these to observations from 
the Hubble Space Telescope. We find the radial distributions that best describe the observed H 2 FUV 
emission arising in primordial disk targets (full dust disk) are demonstrably different than those of 
transition disks (little-to-no warm dust observed). For each best-fit model, we estimate inner and 
outer disk emission boundaries (r in and r out ), describing where the bulk of the observed H 2 emission 
arises in each disk, and we examine correlations between these and several observational disk evolution 
indicators, such as ni 3 _ 3 i, ri ni co, and the mass accretion rate. We find strong, positive correlations 
between the H 2 radial distributions and the slope of the dust SED, implying the behavior of the 
molecular disk atmosphere changes as the inner dust clears in evolving protoplanetary disks. Overall, 
we find that H 2 inner radii are ~4 times larger in transition systems, while the bulk of the H 2 emission 
originates inside the dust gap radius for all transitional sources. 

Subject headings: stars: pre-main sequence - circumstellar matter - protoplanetary disks - molecules - 
line: profiles - ultraviolet: stars 


1. INTRODUCTION 

Protoplanetary disks (PPDs) provi de the r aw materials 
for the formation of stellar systems dBrown et all 120091 : 
IWoitke et al.ll2009l : iDullemond fe Monniedl2010D . Planet 
formation occurs near the midplane of a PPD, w here col¬ 
umn densities and optical de pths are high (jTrilling et al.1 
l2002t lArmitage et al.l 120031 1 , making it difficult to di¬ 
rectly o bserve the material inv olved in the formation 
process ijKominami fe Idal 120021 ). Current understand¬ 
ing of the formation and evolution of planetary sys¬ 
tems in gaseous disks comes from studies of molecular 
content above or near disk midplanes, which place lim¬ 
its on the composition and density distribution of the 
gas and dust cont ent in the inner (r < 10 AU) planet- 


forming regions (Agundez et al.l 

2008): Carr & Naiita 

2001 20llf Salvk et ali 20081 

2011al: IWoitke et al. 

2009; Willacy fc Woods 200l E 

einzeller et al. 20lH: 

Naiita et al. 2011|). “Transition” disks refer to a 


class of PPDs with an optically thick outer zone but 
an inner region significantly depleted of dust grains 
dSato fe N akagawal Il999t I Cal vet et al.ll2002t ISalvk et~ahl 

I2009T ) and are traditionally identified by the deficiency 
in near-infrared (IR) fl ux and steep rise of far-IR flux 
in the obs e rved SED dStrom et all 119891 : I Cal vet et, al.l 
12002L [2005; Es naillat et al.ll2007al) . Several theories ex- 
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ist for how dust gaps are op ened in transition disks , 
including photoevapor a tion dHollenba ch et al. l 11994 
Alexander et alJ 120061 120131 lAlexander fc Armitagel 
20071 IGorti eTffl I2009I). dynamical c l earing by proto- 


planetary systems dCalvet et al.l l2002t iRice et all 12004 
iDodson- Robinson fe Salvkll201lD . and dust grain growth 
ijTanaka et al.ll2005l ). 

Discoveries of significant quantities of gas left within 
the d ust gaps of tran s ition disks (see jNaiita et al.l 
120031 : iRettig et al.l 12004 ISalvk et~al1120071) and sharp 
“walls” betwe e n th e thin and thick dust disk regions 
(jBrown et al.ll2008i) support the possibility of transition 
disks being carve d out b y giant planet formation a nd 
evolu t ion (ISalvk et alJ 20091 : iDodson-Robinson fc Salvkl 
120111 : iDong et al.l 12014 ). The remnant gas disks pro¬ 
vide constraints on the processes that create the fi¬ 
nal structure of planetary systems, such as the transfer 
of gas from the PPD to circumplanetary disks, poten¬ 
tially leading to growth of prot o planets dLubow et al l 
19991 iLubow fc D Angeld 120001 : lAvliffe fc Batel 120101 : 


Beck et al . 2013. Additionally, the molecular atmo¬ 


sphere of transition disks may respond to the dynamical 
perturbations caused by the presence of giant protoplan¬ 
ets and can lead to potentially observable effects, such 
as line asym metries an d d istorti ons in near-IR CO emis¬ 
sion profiles (jRegalv et al.l f2010). The strength of molec¬ 
ular emission originating from the inner radii of PPDs 
is dependent on the gas tempera ture, density, and de¬ 
gree of grain growth (|Salvk et al.ll20lfq) . Molecular line 
surveys therefore provide the opportunity for a broad 
examination of the gas di stributions in circumstellar en¬ 
vironments dBrown et alll20i3 f). 

Molecular hydrogen (H 2 ) has been measured to be 
~ 10 4 times more abundant than any other molecule 
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in the inner disks of PPDs (|France et al.l l2014al) . 
Depending on the density, Ho can surviv e at tem¬ 
peratures up to 5000 K (jWilliamd 12000T ) and self¬ 
shields against UV radiation, mak ing it robust to bo th 
collisional- and photo-di s sociat ion (IBeckwit h et al.lll978t 
IBeckwith fe Zuckermanl Il982t IBeckwith et al.l 11983ft . 
Molecular hydrogen provides a diagnostic for the spa¬ 
tial and struct ural extent of the warm molecular sur¬ 
face of PPDs (lArdila et all 120021 : IHerczeg et al.l 12004 
lYang et all 120111) . While photo-excited H 2 does not 
interact strongly with evolving protoplanets, it traces 
the u nderlying dis t ributi o n of gas at p l anet-f o rming 
radii (lArdila et all l2002t IHerczeg et al.l 12004 120061 : 
IFrance et alJ2012al) . However, H 2 has proven difficult to 
observe in PPDs: cold H 2 (T ~ 10 K) does not radiate 
efficie ntly because it has no permanent dipole dSternberal 
EMIi), so IR ro-vibrational transitions are weak, making 
them difficult to observe from the ground. Therefore, 
studies of molecular material in disks typically rely on 
other tracers available in the near- and mid-IR, such as 
CO and H 2 0, to estimate the molecular disk environ¬ 
ment and mass of the underlying H 2 reservoir in disks. 

The strongest transitions of H 2 are found in the 
FUV (912 - 1700 A), where dipole-allowed electronic 
transitions are primarily photo-excited (“pumped”) by 
Lya photons generated near the protostella r surface 
({France et al.1 l2012bl: iSchindhelm et al.l I2012H) . Warm 
H 2 (T > 1500K) has a significant population in excited 
vibration (v = 1, 2) and rotation quantum states of the 
ground electronic band (X x E+) ( Shul l 1797 8b When a 
Lya photon interacts with a warm H 2 molecule in the 
correct ground-state population [z>,J], the H 2 molecule 
absorbs the photon, exciting it to vibration levels (v' —> 
0-4) of the first electronic band (B 1 E+). Si nce molec ular 
hydrogen has strong (A„; ~ 10 8 s -1 ; see lAbgrall et all 
ll993f ) electronic transitions in the FUV, the excited H 2 
“immediately” decays back to the ground state, emit¬ 
ting a fluorescent photon, observed as an FUV emis¬ 
sion line. The probability for an H 2 excitation-to-ground 
state transition to emit a photon with wavelength A de¬ 
pends on the branching ratio of the allowed transitions to 
the ground electronic state. The brightest H 2 emission 
lines arise from excited states [?/, J'] = [1,4], [1,7], [0,1], 
and [0,2], which have absorption coincidences with Lya 
within 0 and 600 km s” 1 of the Lya line center, large os¬ 
cillator stren£thSj_andrelatiyely_low energy ground-state 
levels (IHerczeg et al.l 120021 1200511 . The set of emission 
lines produced in the [i/, J'] —> [v", J"] decay is refered 
to as a progression. 

Previous work on FUV fluorescent H 2 emission utilized 
basic profile fitting or small-sample parametric sets to 
estimate inner disk diagnostics, such as column density 
and temperature of the r adiating molecular p opulations 
(see IHerczeg et all 12004 IFrance et all 1 2012 al ibi) . In this 
study, we create 2D radiative transfer models of PPD 
atmospheres to reproduce observed FUV H 2 emission 
lines. The models simulate a disk with radial temper¬ 
ature and density distributions, which depend on phys¬ 
ical parameters of the stellar system, such as the disk 
incli nation angle and stellar L ya radiation profile (taken 
from [Schindhc lm et al.ll2012bD . Using the four strongest 
H 2 progressions, we compare radiative transfer emission 
line models to the spectra of 14 CTTSs (8 primordial, 


6 transition disks) observed with the Hubble Space Tele¬ 
scope (HST )/Cosmic Origins Spectrograph (COS) and 
Space Telescope Imaging Spectrograph (STIS). The goal 
of this modeling work is to examine the relationship be¬ 
tween the evolution of warm dust in PPDs and the radial 
distribution of H 2 in the disk atmosphere. We aim to un¬ 
derstand how the spatial distribution of warm H 2 relate 
to the structure of the dust disk and other well-studied 
molecular disk tracers, such as carbon monoxide (CO) 
and water (H 2 0). 

In §2, we present the targets, observations, and se¬ 
lection criteria of H 2 emission features explored in this 
work. In §3, we describe the forward modeling process 
for estimating the warm H 2 disk radiation fields, and in 
§4 we analyze how the best-fit models are determined 
and define metrics used to quantify the evolution of H 2 
radiation for each PPD. In §5, we discuss how the mod¬ 
eled radiation distributions of fluorescing H 2 evolve in 
PPDs, comparing our results with observable warm dust 
disk evolution, mass accretion rates, and additional in¬ 
ner disk molecular tracers. Finally, we summarize how 
the gas disk structure correlates with the dissipation of 
warm dust grains as PPDs evolve to debris disks in §6. 

2. OBSERVATIONS AND H 2 EMISSION LINE SELECTION 


We sample a large collection of HST-COS and HST- 
STIS (for TW Hya) FUV H 2 data to understand the 
relative changes in the radiation distributions of H 2 aris¬ 
ing from the inner regions of primordial and transition 
disks. The observations were obtained through the DAO 
of Tau guest observing program (PID 11616; PI - G. 
Herczeg), the COS Guaranteed Time Observing program 
(PIDs 11533 and 12036; PI - J. Green), and HST Pro¬ 
gram GTO-8041 (PI - J. Linsky). The observations have 


Herczeg et al.l 

20061; Inglcby et al.1 2011;'Yang et al.H2011 

France et al.l 

2012U 2014bl Schindhclm et al. 2012a 


The medium -resolu tion G130M and G160M FUV 
modes of COS dGreen et al.l l2012h were utilized for all 
targets except TW Hya, which was observed with the 
E140M mode (1170 - 1710 A) with the 0.5" x 0.5" 
aperture o f STIS at a resolving power of 25,000 (see 
IHerczeg et al.1 1200611 . The point-source resolution for 
each mode on COS is Av ss 17 km s -1 wi th 7 pixels 
per resolution element (|Osterman et al. 2011) and Av « 
12 km s -1 for the STIS E140M observing mode of TW 
Hya (ILeithered 1200111 . The COS data were smoothed 
by 3 pixels for analysis. The one-dimensional spectra of 
COS were produced using the CALCOS COS calibration 
pipeline, which were aligned and coadded u sing a custom 
software procedure (iDanforth et al.l 120101) . The STIS 
data wer e reduce d using the CALSTIS STScI reduction 
pipeline (ILindlerl 119991) . with calibration lamp spectra 
obtained during observations to assign wavelength solu¬ 
tions. An example of the continuous far-UV spectrum of 
V4046 Sgr is shown in Figure [l] 

Stellar properties, such as mass, accretion rate, and 
inclination angle are used to constrain the underly¬ 
ing model framework. All disk inclination angles 
have been estim ated from sub -mm/IR in terferometric 
studies (see iSimon et al.1 20001: Johns-Krull fe Valent 
20011: lAndrews fc Williams 20071: Espaillat et all 12007: 


Andrews et al.l 120111 : iRosenfeld et al.1 I2012D . Stel- 
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Table 1 

Stellar Parameters 


Target 

Spect. 

M* 

d. A-y 

id 

Age 

vsini 

ref. a 


Type 

(M 0 ) 

(pc) 

n 

(Myr) 

(km s _1 ) 



AA Tau 

K7 

0.8 

140 

0.5 

75 

6.4 

± 

0.2 

11.4 

2,8,9,11,12,15,17,25 

BP Tau 

K7 

0.73 

140 

0.5 

30 

5.9 

± 

0.3 

7.8 

4,8,9,12,15,17,18,30 

CS Cha 

K6 

1.05 

160 

0.8 

60 

6.4 

± 

0.1 


5,6,15,19,22 

DF Tau A 

M2 

0.19 

140 

0.6 

85 

6.3 

± 

0.5 

16.1 

9,11,12,16,17,18 

DM Tau 

Ml.5 

0.5 

140 

0.0 

35 

6.6 

± 

0.2 

10.0 

3,11,12,15,17,25 

GM Aur 

K5.5 

1.20 

140 

0.1 

55 

6.9 

± 

0.2 

12.4 

3,8,9,11,12,15,17,25 

HN Tau A 

K5 

0.85 

140 

0.5 

40 

1.9 

± 

0.9 

52.8 

7,9,11,18,23 

LkCal5 

K3 

0.85 

140 

0.6 

49 

6.4 

± 

0.3 

12.5 

3,8,10,12,15,17,18 

RECX 11 

K4 

0.80 

97 

0.0 

70 

4.0 

± 

1.5 


14,15,20,21 

RECX 15 

M2 

0.40 

97 

0.0 

60 

6.0 

± 

1.0 


15,20,21,32 

SU Aur 

G1 

2.30 

140 

0.9 

62 

2.5 

± 

0.9 

65.0 

1,4,9,11,18 

TW Hya 

K6 

0.60 

54 

0.0 

4 

10.0 

± 

6.0 

6.0 

3,13,16,24,28,29,31 

UX Tau A 

K2 

1.30 

140 

0.2 

35 

6.1 

± 

0.3 

25.4 

3,8,11,17,19 

V4046 Sgr 

K5 

1.75 

83 

0.0 

34 

6.9 

± 

0.1 

14.2(+13.7) 

25,27,28,29 


(4) Bouvierl ( 

1990); (5) lEs 

paillat et al. (2007a); (6) Esi 

:>aillat et al. ( 20111); (7) France et al. 

(20121 

); (8) 

Furlan et al. 


2011 

1989 

); (9) Gullbring et al. 

199g); (10) Hartmann et al. 

(1987); 

(11) H 

Lartmann &; Stauffer 

); (12) Hartmann et al. 

(1998): (13) llerczee L Hil 

enbrand 


I Lawson etlr l. (2001); (21) Lawson ct al. 1120041b 
iPo ntoppidan ct aTI |(200S|^ 


_ (22) lTmlnna.nl 1120041 b (231 IMc.Tunki'n et all 

(25) [t juast et al] <[200Ch: (261 IRicci et HI IMinh : 


(f2(TTab (241 _ ._ _ _ 

(27) Rodriguez ct al. (2010); (28) Rnscnfcld ct al. (2012); (29) Roscnfcld ct al. (2013); (30) 
ISimon et**alTl2'oO0l b (311 IWebb et alj d 199911 : (321 IWoitke et al] (120131 1. 


lar masses and extinction estimates were derived 
from pre-main seque nce stellar evolutionary tracks 
(iHartmann et al.lll998l) . Mass accretion rates were es- 
timated from meas urements of the accretion luminosity 
(llnglebv et al.ll2013h . Refer to Table Q] for lists of all the 
relevant stellar parameters, with references therein. 

The HST-COS FUV spectra of all CTTSs from 1300 
- 1600 A reveal a suite of H 2 fluorescence features linked 
to Lya-pumping. We chose to use the strongest transi¬ 
tions from the electronically-excited progressions \v ', J'] 
= [0,1], [0,2], [1,4], and [1,7] for the purposes of study¬ 
ing the underlying general characteristics of the bulk gas 
disk. We sample 3 emission features from each progres¬ 
sion. This gives us access to strong, non-blended emis¬ 
sion lines that are well-defined from the FUV continuum, 
while balancing the CPU time required for detailed line 
profile analysis. We selected H 2 emission features by lo¬ 
cating the strongest tra nsitions for each progression, out¬ 
lined bv iHerczeg et all (|2002[ ). See Table [I] for the full 
outline of transitions chosen. 

3. MODELING ANALYSIS 

We create models of warm H 2 in PPDs to constrain 
the radial distribution of fluorescent H 2 emission in disk 
atmospheres. Our aim is to understand the relative 
changes in the H 2 distributions as we observe various 
stages of dust disk evolution. The fluorescent emission 
line shape and intensity depend on the physical condi¬ 
tions of the gas, while the observed line width depends 
predominantly on the disk inclination. We construct a 
physical model of the disk str ucture, motivated by t he 
disk modeling analysis done bv iRosenfeld et all (|2012f ) . 

The models make several basic assumptions on the 
disk properties: (a) the disk material orbits in Keple- 
rian rotation around a central point mass, representing 
the stellar mass; (b) the H 2 fluorescence occurs in a shal¬ 
low, warm layer on the disk surface; and (c) the level 
populations of warm H 2 that absorb the incident stel¬ 


lar Lya radiation field are in local thermodynamic equi¬ 
librium (LTE). (a) implies that the gas disk mass is a 
small fraction of the stellar mass (Af^/Af* -C 1). Several 
studies have shown that the disk mass to stellar mass 
ratio (Md/M+) < 1% , making this assumption plausible 
([Andrews et al.l I2013D . In the case of a binary system 
(i.e., V4046 Sgr), both stellar mas ses are represe n ted as 
one central mass point. For (b), IHerczeg et al.1 (|2004f ) 
find that the warm H 2 disk layer interacting with the 
stellar Lya to produced the observed fluorescence lines 
corresponds to mass column density of ~ 10 -5 g cm -2 , 
which is a much smal ler mass column densit y predicted 
to be within 1 AU by iD’Alessio et al.1 (|l999t ). This sug¬ 
gests that the Lya-pumped fluorescent emission origi¬ 
nates from a tenuous layer of warm H 2 on the disk surface 
and supports a purely radial thermal distribution T(r). 
For (c), combination of collisional excitation and radia¬ 
tive de-excitation is assumed to be in equilibrium to keep 
the H 2 gas n ear the disk surface at wa rm temperatures 
(T > 1000 K dNomura et al.ll2005t 120071) . Previous stud¬ 
ies of FUV H 2 e mission have argued both for and agains t 
this assumption (lArdila et al.l[2002l IHerczeg et al .11200(1) . 
LTE conditions keep the assumed parameters straightfor¬ 
ward and allow us to model the H 2 ground-state popu¬ 
lations as a “snapshot” of the disk atmosphere as it was 
observed. 

The warm H 2 atmosphere is described by the surface 
density and temperature distribution of gas, which char¬ 
acterizes how much of the warm H 2 is populating excited 
ground-states [v,J\. We reference these physical quanti¬ 
ties in cylindrical coordinate positions in the disk (r,<j>,z). 
If we consider that a parcel of warm H 2 gas on the disk 
surface is characterized by its radial position, vertical 
height from the disk midplane, and velocity distribution 
(r, z, v<j,(r)), the velocity of the gas parcel, v^r), is de- 
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1300 1350 1400 1450 o 1500 1550 1600 

Wavelength (A) 

Figure 1. The HST -COS FUV spectrum of transition disk target V4046 Sgr, ranging from 1300 - 1600 A. Most of the narrow emission 
features are Lya pumped fluorescent H 2 lines. Many of the strongest H 2 emission features fluoresce from Lya-pumped excited state 
progressions [1,4], [1,7], [0,1], or [0,2]. Blue diamonds mark the H 2 fluorescent emission lines studied in this work. 

scribed by Keplerian rotation in </> only: of the disk atmosphere is approximated as a power-law 

_ function: 

M r ) =Vk = =v z =0, ( 1 ) T(r) = T\ AU ? . ( 4 ) 


where G is the gravitational constant and M* is the cen¬ 
tral stellar mass. The mass density at the warm H 2 disk 
surface is a function of the radial and vertical height in 
the disk, 


p{r, z) 


£Q) 

V^H p XP 2 UJ 


( 2 ) 


where £ (r) is the radial surface density distribution of 
H 2 , and H p is the pressure scale height as a function of 
radius, defined as: 


r c s lkT(r) r^~ 

p ~ fl ~ Y fj,m H ' GM * ’ 


(3) 


where c s is the sound speed, is the angular velocity of 
the gas, k is the Boltzmann constant, T(r) is the radial 
temperature profile of the warm H 2 disk atmosphere, /i 
is the “mean molecular weight” of the gas, and to h is the 
mass of a hydrogen atom. The temperature distribution 


where T\au is the temperature of the warm H 2 at r = 1 
AU, and q is the temperature gradient. 

We assume a radial surface density for a static accre¬ 
tion disk , represented by a pow er-law viscosity profile 
(see iLvnden-Bell fe Pringleifl974ll , 


£(r) = £ c 




(5) 


where 7 is the density gradient, r c is the characteris¬ 
tic radius of the gas in the disk, and £ c is a normal¬ 
ization factor for the surface density distribution, de¬ 
pendent on the total H 2 mass contributing to the emis¬ 
sion lines simulated by these models. The characteristic 
radius describes the transition from a power-law domi¬ 
nated density distribution to an expo nent ially-dominated 
density fall-off in the disk (|Lvnden-Bell fe Pringlel[l974 
lHartmann et alJll998l) . It is important to note that £(r) 
contains a normalization factor (£ c ), which normalizes 
to the disk midplace density. Our models only attempt 
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Table 2 

Selected H 2 Emission Lines & Properties of H 2 Pumping Transitions 


A lab 

Progression 

Line ID a 

A pump 

Vtrans 

A u i c 

/ d 

(A) 




(A) 

(km S -1 ) 

(10 s s 1 ) 

(10- 3 ) 

1442.87 

[1,7] 

(1 

- 6)R(6) 

1215.726 

14 

0.9 

34.8 

1467.08 


0 

- 6)P(8) 



1.3 


1500.45 


(1 

- 7)R(6) 



1.7 


1524.65 


(1 

- 7)P(8) 



1.9 


1556.87 


(1 

- 8)R(6) 



1.3 


1580.67 


(1 

— 8)P(8) 



1.1 


1431.01 

IMI 

(1 

- 6)P(3) 

1216.070 

99 

1.0 

28.9 

1446.12 

(1 

— 6)P(5) 



1.4 


1489.57 


(1 

- 7)P(3) 



1.6 


1504.76 


(1 

- 7)P(5) 



2.0 


1547.34 


(1 

- 8)P(3) 



1.1 


1338.56 

[00] 

(0 

- 4)P(2) 

1217.205 

379 

3.1 

44.0 

1398.95 

(0 

- 5)P(2) 



2.6 


1460.17 


(0 

- 6)P(2) 



1.5 


1521.59 


(0 

- 2)P(2) 



0.6 


1342.26 

[0O] 

(0 

-4)P(3) 

1217.643 

487 

2.8 

28.9 

1393.96 


(0 

- 5)P(1) 



1.6 


1402.65 


(0 

- 5)P(3) 



2.3 


1463.83 


(0 

- 6)P(3) 



1.4 


1525.15 


(0 

- 7)P(3) 



0.5 



a Transitions are from the Lyman-excited to ground electronic states of the H 2 band 
system, B 

k Velocity from line center of the pumping transition of Lya. 

c Einstein coefficient, describing the spontaneous decay rate from the electronically- 
excited Lyman band, taken from (Tlmral I et all (11993 !■ 

^ Oscillator strengths from lAbgrall et all (1 199311 - 


to describe the behavior of the disk atmosphere, where 
the warm, tenuous H 2 resides. As a consequence, the 
functionality of E(r) serves as a structural layout of the 
radial H 2 disk atmosphere. Since we normalize X(r) with 
a factor describing the disk midplane density, the solu¬ 
tions of E(7') describe the radial distributions of warm 
H 2 , but the resulting H 2 mass estimates are not mean¬ 
ingful. 

The level populations of warm, ground state H 2 con¬ 
tributing to the emission line are assumed to be in LTE 
and are determined using the Boltzmann equation, 


n[ Vl j]{r,z) 


p(r,z) X H2 
pm H 


S[-U, J] 

Z [v,J] ( T ) 


x exp 


U T{r)J' 


( 6 ) 


where X# 2 is the fraction of the total H 2 gas mass con¬ 
tributing to the fluorescence observed in the FUV, g[u, j] 
is the statistical weight of the level population, Zu, n(T) 
is the partition function describing the likelihood that 
the warm H 2 is in state and E^ j] is the energy of 

warm H 2 in ground state [u, J], 

The radial distribution of molecular hydrogen has two 
normalization factors (Xh 2 and E c ) that are not in¬ 
dependent of disk conditions and are defined by their 
product in n\ v ,j] (r, z). The product of these factors de¬ 
scribes the total mass of warm H 2 available for photo¬ 
excitation to state [v',J'] ( Mh 2 ), which is obtained 
by integrating the distribution over ( [r,<f>,z ): Mh 2 = 
X* 2 £ c (27rr c 2 )/(2-7). 

The radiative transfer calculation required to re¬ 
produce the observed fluorescent H 2 emission happens 
in two steps: 1) the warm H 2 in ground state population 
[.X : v,J] is pumped into a rovibrational level [B : v',J'] of 



0 2 4 6 8 10 


r (AU) 

Figure 2. We provide a graphical representation of the H 2 disk 
atmosphere model. The disk contours represent the warm H 2 opti¬ 
cal depth (r\(r, z)) to stellar Lya radiation being pumped to state 
= [1,4]. The dashed line marks off the approximate loca¬ 
tion of Ty ~ 1, which is where the H 2 disk atmosphere becomes 
optically thick to the penetrating Lya photons. The stellar Lya 
radiation (purple arrow) is absorbed by the by the warm H 2 , which 
is excited to state [v',J'] and emits a photon (A# 2 ; red arrow) to 
decay back to ground state [v" ,J"]. 

the excited electronic (Lyman band) state by the absorp¬ 
tion of an incident stellar Lya with wavelength A Lya , and 
2 ) the excited H 2 molecule decays back to some ground 
electronic state [X : v" ,J"], emitting a FUV photon 
with wavelength A h 2 - Molecular hydrogen has an ab¬ 
sorption cross section ( cth 2 ) defined by the area around 
the molecule that can intersect an incoming photon with 
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Table 3 

Parametric Values Explored in Modeling Framework 


Table 4 

Minimum x 2 Statistics for Each Progression Fit 


Parameter 

Values 

Units 

z/r 

(2, 3, 4, 5, 6 , 7) x H p 


7 

0.0, 0.25, 0.5, 0.75, 1.25, 1.5, 1.75, 1.99 


q 

-1.0, -0.5, -0.25, -0.1, -0.05, 

0.0, +0.05, +0.1, +0.25, +0.5 


Tl AV 

500, 1000, 1500, 2000, 2500, 

3000, 3500, 4000, 4500, 5000 

K 

r c 

0.1, 0.5, 1.0, 3.0, 5.0, 7.5, 10.0, 20.0 

AU 

M H2 

5xlO~ 10 , 10“ 10 , 5xl0 -11 , 

10“ u , 5xl0 -12 , 10“ 12 , 
5xl0 -13 , 10 —13 , 5xl0 -14 , 10 -14 

M 0 


Values were chosen to reproduce the desired H 2 features 
dHerczeg et al.| I2004J ; IFrance et al.| |2Q12bf) . The only parameters 
without aforementioned constraints were z/r, 7 , q and r c because 
literature estimates of these values were not known. 7 and q 
were constrained by the power-law functionality role they play in 
the models, and r c was estimated around <r[i 5 7 ]> calculated by 
IFrance et al.| d2012bD . 

the appropriate energy for photo-excitation: 


0/2 


^Lya S[B:v',J'] . 

s - A lui 

8/rc g [X:v ,j] 


(7) 


where A Lya is the rest frame wavelength of the stellar 
Lya line profile needed to excite the warm H 2 in ground 
state [X : v, J] up to energy level [B : v ', J'], and Ai u is 
the probability that H 2 in population [A' : v , J] will be 
“pumped” to electronic state [B : Note that, for 

the remainder of this paper, we will omit the ground state 
H 2 ( X ) and excited state ( B ) level branch denominations 
from the vibration and rotation state discussion. 

Assuming an absorption coefficient K\ ( r , z) = 
(7 H 2 n \v,j\ ( r > ~)> the optical depth of H 2 in ground state 
[v,J] is described as: 


z Hp 

t\ ( r , z) = ^2 ZKx (r, z ). (8) 

Z 


For every vertical and radial position in the disk atmo¬ 
sphere that we sample T\(r, z), we calculate the amount 
of the Lya radiation that will be available for absorption 
by the warm H 2 . To correct for line absorption overlap of 
shared Lya photons, we adopt an effective optical depth 
t' x (r, z) dLiu fe Dalearnol 119961 : iWolven et al.l I1997D , de¬ 
fined as 


T \ (A z) = t\ ( r , z) 


t\ (r, z) 


^11 (T(r), N(r, z)) ’ 


(9) 


which corrects for the absorption, scattering, and shield¬ 
ing of Lya photons. Figure [2] shows a schematic of 
T\(r, z) for [v', J'} = [1,4] and outlines the radiative trans¬ 
fer process in the disk. 

We model the emission line flux of each A h 2 produced 
from the cascade of transitions from energy level [v', J'} 
as: 


Fx H2 


r]S\ {r, z) B mn 



( 10 ) 


Target 

[ 0 , 1 ] 

Progression 

[ 0 , 2 ] 

W,J'} 

[1,4] 

[1,7] 

AA Tau (2011) 

5.37 

6.48 

5.52 

1.25 

AA Tau (2013) 

1.78 

5.29 

4.24 

1.62 

BP Tau 

2.82 

51.75 

5.28 

2.97 

CS Cha 

4.56 

5.14 

4.19 

2.62 

DF Tau A 

2.69 

13.30 

7.21 

7.37 

DM Tau 

6.12 

19.55 

7.68 

37.95 

GM Aur 

3.84 

6.72 

1.47 

1.74 

HN Tau A 

41.71 

63.16 

13.52 

35.81 

LkCal5 

111.03 

103.30 

14.14 

151.65 

RECX 11 

2.40 

9.48 

1.09 

0.93 

RECX 15 

42.45 

90.01 

13.98 

63.32 

SU Aur 

25.73 

39.31 

13.24 

21.07 

TW Hya 

2.64 

3.29 

3.63 

2.15 

UX Tau A 

104.69 

124.23 

13.14 

123.16 

V4046 Sgr 

12.82 

13.09 

5.93 

2.86 


All model-to-data reduced-x 2 statistics for simultaneous 
emission line fitting, transitioning from excited state [v r ,J'\. 
All x 2 statistics are calculated between v 0 b s = [-250, 250] 
km/s. The largest source of errors in the x 2 statistics come 
from the linear estimation of the FUV background continuum 
beneath the emission line. Because the models do not attempt 
to find the background continuum levels beneath each emis¬ 
sion line, extraction of the FUV continuum had to be done 
manually. Targets with lower signal-to-noise have more uncer¬ 
tainty in the continuum flux, so the x 2 statistics become large 
as the errors in the continuum dominate the fitting. Only the 
[1,4] progression show decent fits for all targets (with x 2 < 
15), so we focus on the relative results of the [1,4] progression 
emission lines for the remainder of the Discussion section. 



Figure 3. An example of a modeled emission line fit over a HST- 
COS emission line. The target is DM Tau, one of the transitional 
disk targets, for emission line A1489.57 A, which fluoresces from 
the Lyman band energy level [1,4]. The black line represents the 
observed H 2 fluorescent emission feature, which includes error bars 
every 5 bins. Each emission line observed has an intrinsic back - 
ground continuum from the stellar source (see lFrance et al.l2014bTh 
so this continuum was subtracted from the line before model com¬ 
parisons were made to the observations. The red line is the modeled 
emission of 1489.57A from the DM Tau disk model. The blue line 
is the convolution of the modeled emission line with the COS LSF. 
This procedure was applied to all modeled emission lines for all 
targets when comparing the modeled data with FUV observations. 
The reduced-x 2 was calculated after the model emission lines were 
convolved with the COS LSF. 

where rj represents the coverage fractio n of H 2 in the 
Lya radiation field dHerczeg et al.1 12004H . B mn is the 
branching ratio describing the fraction of H 2 decaying 
via a given transition to ground state \u" ,J"\ over the 
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Figure 4. The resulting model and data fits of the minimum simultaneous progression y 2 statistic for GM Aur. Each column represents 
transitions from a common excited energy level \v',J']. From left to right: the left column - \v', J'\ = [1,7]; the middle-left column - [1,4]-; 
the middle-right column - [0,1]; the right column - [0,2]. All reduced-x 2 values for each progression were calculated by simultaneously fitting 
each observed emission line profile to those estimated in a given model parameter set. The minimum reduced-x 2 for each progression is 
assumed to best represent the H 2 fluorescence distribution. The x 2 shown in the top left of each emission line box represents the fitting 
of all emission lines from a gi ven progressio n with one set of model parameters. Most of the observed emission lines for all targets have 
single-peaked line profiles (see I Brown et aLll2013i 'l. but all the best-fit modeled emission lines show a “double horned” profile signature to 
Keplerian motions. IPontoppidan et al.l 112011 :'! points out that the single-peaked profile behavior is suggestive of a low velocity (< 3 km/s) 
molecular wind located within a few AU of the central star and are typically modeled with an azimuthal velocity vector that is slow relative 
to Keplerian motion (also see lBast et al.|[2Q] I ). Since our models leave out the azimuthal velocity component of the H 2 disk gas motions, 
it is expected that our resulting emission line profiles do not reproduce the line cores of the fluorescent features. 


whole suite of transitions available from the progression, 
and the source function (S\(r,z)) is defined as the Lya 
emission line flux with wavelength \L ya , F^ ya (r,z). 

We calculate how Fi, ya (r,z) changes as a function of 
radial position in the disk. Assuming that the accretion¬ 
generated Lya flux originates at the stellar surface, we 
express the ratio of the original Fj AjQi , to the flux the 
warm H 2 disk atmosphere receives at r, 


F, 


Lya 


= F 


Ri 


*’ L V a r 2 


( 11 ) 


To correctly incorporate the Lya radiation field, 
we use reconstru c ted st ellar L va profiles cr e ated b y 
ISchindhelm et al.l (I2012bl) and iFrance et al.1 (|2014bl) . 
which describe the stellar-Lya flux seen by the disk sur¬ 
face of each target. After calculating the FUV H 2 flu¬ 
orescence flux at each disk grid point in our model, we 
radiate the H 2 emission isotropically, some fraction of 


which is intercepted by the observer. We calculate the 
distance of each gas parcel radiating in the disk from the 
observer s(r,z), based on radial and angular positions 
of the disk gas parcel, distance to the target, and disk 
inclination angle. The final modeled emission line flux 
produced for a fluorescence transition of H 2 is expressed 
as: 


F\h 2 



(d cos i d isk) 2 \ 
s{r,z) 2 J X 



( 12 ) 


Using a total of 6 parameters to represent the physical 
conditions of the warm, ground-state H 2 populations in 
the disk atmosphere (z/r, ”/,q,TiAUii"char, M# 2 ), Equa¬ 
tion m characterizes the resulting emission line profiles 
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from H 2 radiating from the disk. All free parameters were 
allowed to vary over a rough grid of controlled values to 
create a data cube representing the density distributions, 
temperature profiles, and radial radiation fields of inner 
disk H 2 around a given stellar target; see Table [Q for the 
full list of parameters explored in this study. The result¬ 
ing models simulate the emission profiles produced for a 
given fluorescence transition A# 2 , with emission flux as 
a function of orbital velocity. The radial velocity com¬ 
ponent of the emission line is determined by v^r) of the 
emitting gas at a given radius in the disk, projected into 
the sight line of the observer. This model framework was 
used to describe the observed velocity field of single and 
binary systems, both close-in and extended. We caution 
the reader regarding the results of the close-in binary sys¬ 
tems (e.g. V4046 Sgr), as the binary potential affects the 
inner disk velocity-radial relationship differently than a 
point mass. Therefore, the innermost H 2 modeled for 
these close-in binary systems may not be accurate, but 
the outer disk emission distributions will remain unaf¬ 
fected. 

Synthesized spectra of each H 2 emission line are 
compared to HST observations. Each model is con¬ 
volved with eithe r the HST-COS line spread function 
(LSF) (1 K i~ is sll 2 (iTTI) or a normalized Gaussian distribution 
with FWHM characterized by the STIS E140M mode 
spectral resolving p ower (R/^25,000 for TW Hya; see 
IHerczee et al.1 [20061 for more information) prior to com¬ 
parison with the observed emission line profiles. The 
FUV continuum level is estimated around each emission 
feature with a linear fit to the HST- COS data, which is 
subtracted from the observations before model-to-data 
comparisons are made. An example of an H 2 emission 
line, with native and convolved models laid over the 
HST- COS observed emission line, is shown in Figure [3] 

4. ANALYSIS 

The goal of the model-to-data comparison is to find 
the combination of model parameters that best repro¬ 
duce the observed fluorescent emission line profiles that 
cascade from the same excited state [v',J']. A reduced-y 2 
statistic is computed when comparing the observed FUV 
H 2 emission features to the entire data cube of models 
created for a target. We analyze the reduced-y 2 statistic 
data cube for three cases when comparing the modeled 
emission lines to the observations: (1) fitting individ¬ 
ual emission lines; (2) simultaneously fitting all H 2 emis¬ 
sion lines fluorescing from the same excited energy level 
[ 1 /, J']\ (3) fitting only the red wings of the emission lines. 
The first point was used to set the initial range of temper¬ 
ature and density model parameters of warm H 2 in each 
disk surface. The third was explored to mitigate the po¬ 
tential influence of a warm molecular wind component 
that was unresolved at the spectral resolving power of 
HST-COS. The results of (3) proved inconclusive, which 
found no significant differences between the red and blue 
wing line shapes, suggesting that the models are not sen¬ 
sitive to an unresolved warm H 2 disk wind. We focus on 
the results of (2), which best describe the generalized be¬ 
havior of the warm H 2 disk atmosphere populations. We 
simultaneously fit 3 observed fluorescent H 2 transitions 
for each progression as the most representative of the H 2 
radiation distributions in each PPD. 

Table □ shows the minimum reduced-y 2 statistics for 



Figure 5. Using the best-fit progression model for GM Aur, we 
use Equation 12 (integrated over B mn , which represents the total 
H 2 flux produced from each progression) to reproduce the observed 
spectrum. Each progression peaks at different radii, but the overall 
shape and radial extent of the distributions indicate that the bulk 
of the radiation for all progressions originates within the same disk 
annuli. 

all targets when simultaneously fitting the 3 progression 
emission lines from excited state \v' Not all mini¬ 
mum reduced-y 2 simultaneous progression fits for [0,1] 
and [0,2] were “good”, however (i.e., some sources dis¬ 
played reduced-y 2 > 25). Many of the strongest lines 
from [0,1] and [0,2] share similar A# 2 , which makes com¬ 
plex line profiles that depend on the shape of the stellar- 
Lyct profile illuminating the warm H 2 disk populations to 
these excited states. The [1,7] and [1,4] progressions are 
more reliable tracers of the warm H 2 disk atmosphere, 
and the brightest emission lines in our survey cascade 
from the [1,4] progression. For the vast majority of the 
targets, the largest stellar Lya fluxes pump the warm 
H 2 disk populations to the [1,4] energy level. This makes 
the line profile flux fitting more accurate for the [1,4] 
progression, providing the overall best model fits to the 
observe FUV emission. 

We will focus our discussion around the inner disk 
diagnostics of the best-fit [1,4] progression for all tar¬ 
gets. This progression has good reduced-y 2 fits (< 15) 
and by-eye model-to-data comparisons for every target 
in our survey. Figure [4] shows an example of minimum 
reduced-y 2 modeled progression lines to those observed 
with HST-COS for GM Aur. Figure [5] presents the re¬ 
sulting radial radiation distribution for each best-fit pro¬ 
gression for GM Aur. While each progression peaks at 
somewhat different radii, the majority of the radiation 
distributions originate within similar annuli of the disk. 
This behavior is typical for all PPD targets that have 
good minimum reduced-y 2 fits for all or most progres¬ 
sions. 

4.1. Uncertainty Estimation and Parameter 
Degeneracies 

Errors in each best-fit parameter per progression are 
determined after marginalizing the minimum reduced- 
y 2 parameter fits over all free parameters. Uncertainties 
are measured as the range of values that encompass 68% 
of the distribution area, representing the 1 -a uncertain¬ 
ties for a Gaussian distribution. The modeled parameter 
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M H2 (M 50| ) r char|AUI q T , AU (K) y z/r 


Figure 6. Marginalized distributions of the [1,4] progression reduced-x 2 fits for RECX-11. The uncertainties in the best-fit model 
parameters are measured as the range of values that encompass 68 % of the distribution area and are highlighted in red contour outlines. 
The posterior marginalized distributions for each parameter against all other model parameters are shown as the 2D plots at the top of 
each column. A Gaussian distribution was estimated for each posterior distribution, and final uncertainty estimates for each model was 
calculated as the FWHM of the distribution. These errors were later used to estimate errors in radial radiation distribution boundaries. 


space was crudely varied over a large range of values for 
each free variable, so a Gaussian distribution was fit over 
each marginalized best-fit parameter uncertainty space, 
and the FWHM of each Gaussian fit was calculated as 
the uncertainty in each model parameter. 

Figure [6] displays the reduced-^; 2 marginalized para¬ 
metric space for each variable in our modeling frame¬ 
work, with filled contours representing the 2-a uncer¬ 
tainty in the parameter space. Since each parametrized 
uncertainty is taken within the 1-cr error contours of each 
marginalized distribution, the uncertainties outlined in 
red represent the 1-a errors in the model parameters. 

There are noticeable degeneracies amongst several of 
the parameters; for example, the total mass of emitting 
H 2 and vertical position of the disk atmosphere ( Mh 2 , 
z/r) show a trend that requires more mass contributing 
to the emission lines as the disk height above the disk 
mid plane increases. This trend makes sense - to pro¬ 
duce the same amount of flux in the modeled emission 
lines, the total mass of H 2 contributing to the emission 
must increase as the density of H 2 decreases with verti¬ 
cal disk height above the mid plane. The optical depth 
of the disk atmosphere must remain the same to output 
the same observed emission line flux, and this relation¬ 
ship between the free parameters maintains the required 


optical depth. What is important to note is that the 
models produced are used as a means to describing the 
H 2 emission flux arising from the inner disk atmosphere. 
Despite the degeneracies in several parameter pairings 
relating to the total flux, the radiation distribution of 
H 2 emission is unaffected by these degeneracies. 

We note that our choice in using the reconstructed stel¬ 
lar LyA flux incident on the disk from Schind helm et all 
( 201213 ) m ay exacerbate degeneracies in the disk param¬ 
eters. The iSchindhelm et al.l (|2012 b) reconstructed Lya 
profiles rely on the same H 2 emission features explored in 
this study, but we remind the reader that the stellar Lya 
flux incident on the H 2 disk scales with the re-emitted H 2 
flux (see Ean fl2|l and has no effect on the modeled distri¬ 
bution of H 2 flux in each disk. The disk parameters may 
respond to an inaccurate Lya flux, but the degeneracies 
in the disk parameters (for example, the response of Mh 2 
and z/r to the total H 2 flux) would scale to best describe 
the H 2 radiation that recreate the observed emission pro¬ 
files. Therefore, the reconstructed Lya profiles will not 
change the radial behavior of the best-fit H 2 flux models. 

4.2. The Radial Extent of H 2 Emission 

Figure [5] presents an example of the radiation distribu¬ 
tions of H 2 fluorescence flux produced in the disk for each 
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Figure 7. The inner and outer radial boundaries which define 
where 90% of the total radiation is arising from the disk. The black 
line represents the normalized radial distribution of the [1,4] pro¬ 
gression emission to the total amount of flux produced by the [1,4] 
progression for SU Aur. The green vertical lines show the radial 
boundaries that encapsulate 90% of the total emission. For r j n , we 
start at the outermost radius and integrate inward to smaller radii 
the disk until 95% of the total [1,4] progression flux is accounted 
for. Likewise, r ou t is defined by starting at the innermost radius 
defined in our models and integrating the progression flux out until 
95% of the total emission flux is accounted for. The resulting an¬ 
nulus between r i n and r out represents the ring of disk the majority 
of the observed FUV H 2 fluorescent emission originates. 


progression explored in this study. We focus our analysis 
on the [1,4] radiation distributions for all targets in our 
survey to define inner and outer radial disk boundaries, 
which describe where the bulk (90%) of the emitting H 2 
atmosphere resides. We define the 90% emitting region 
as follows: 


F, 


H 2 ,obs 


F(H 2 ,r) 
F to t(H 2 ) 
F(H 2 ,r ) 
F tot (H 2 ) 


< 0.95 for r > r in 

< 0.95 for r < r out 


We use ri n and r ollt to evaluate the evolutionary be¬ 
havior of the H 2 radiation. Figure [7] presents a schematic 
of how the inner and outer radial boundaries encapsulate 
90% of the total H 2 flux produced in the disk atmosphere. 
We analyze potential evolutionary characteristics of the 
molecular disk atmosphere by comparing the FUV H 2 
radiation distributions to other dust and molecular disk 
observables. 


4.3. Case Study: Model Robustness using AA Tau 

We explore how robust our modeling framework is at 
identifying where the fluorescing H 2 resides in PPDs. We 
compare two epochs of HST-COS data on AA Tau (2011 
and 2013), where the 2013 observations occur during a 
“dimming” event from X-ray to near-IR wavelen g ths. 
Based on the duration of the dimming, iBouvier et al.1 
(j2013lf suggest an obscuration at r > 8 AU; this hypoth¬ 
esis is strengthened by the gas-to-dust ratio (N#/A t ,) 
of th e absorber and the evo lution of the FUV H 2 emis¬ 
sion (]Schneider et al.l [2015b We utilize the line pro¬ 
file changes between AA Tau HST-COS FUV observing 
epochs to determine how those changes relate to radial 
H 2 radiation distributions in the disk. 

There are noticeable differences between the observed 
FUV H 2 emission line profiles of the 2011 and 2013 AA 


Figure 8. We present a comparison of the [1,4] progression ob¬ 
served with HST-COS for AA Tau during the 2011 and 2013 
epochs. On the left: The column under the 2011 label represent the 
3 observed fluorescent emission line profiles cascading from the [1,4] 
excited state, with overlaid best-fit modeled emission lines in blue. 
The column to the right, labeled 2013, shows the observed [1,4] 
fluorescent emission lines, with modeled emission lines overlaid in 
orange. The 2013 observed e mission line profiles ap pear narrower 
than their 2011 counterparts dSchneider et al1l2015l) . On the right: 
the comparison of the total [1,4] progression flux radiating from the 
disk of AA Tau in the 2011 and 2013 observations. The 2013 mod¬ 
els predict that the observed H 2 fluorescence emission originates 
from further out in the disk (the peak of the radiation located at 
Tpeak = 2.50 AU) than the 2011 radiation distribution (r pea fc = 0.75 
AU), a consequence of the inner disk “sh adowing” produced b y the 
extra absorber on the AA Tau sightline (IBouvier et al.|[2013D . 

Tau epochs. The 2013 emission lines are narrower with 
slightly larger peak fl uxes than the same H 2 emission 
lines observed in 2011 ((Schneider et al.ll2015H . This sug¬ 
gests that less flux is contributed from the innermost 
disk. The modeling results for the [1,4] progression are 
shown in Figure [H] Each A A Tau epoch was modeled 
independently, and the models reproduce the same rest 
wavelength emission lines. Figure [ 8 ] also shows the radi¬ 
ation distributions of [1,4] fluorescence for each epoch in 
the AA Tau disk. The 2011 emission includes a large con¬ 
tribution from material inside 1 AU (r^^on = 0.08±0.01 
AU; T pe ak, 2 oil = 0.75 AU; r OMtj 2on = 4.17T2.04 AU), 
while the 2013 [1,4] emission “appears” to have shifted 
outward in the disk (^,^2013 = 0.15±0.02 AU; ^^2013 
= 2.50 AU; r o „ tj 20 i 3 = 7.59±2.75 AU). Our models indi¬ 
cate that the inner radius of detectable H 2 fluorescence 
from the [1,4] progression has moved outward radially in 
the disk as the “extra abs orber” moved in t o our field of 
view in the AA Tau disk. iSchneider et al.l (I2015H . using 
an independent modeling technique to estimate the ra¬ 
dial origins of H 2 fluorescence in the AA Tau disk, come 
to a similar conclusion: the observed 2013 H 2 emission 
withi n ~ 1 AU is reduced c ompared to 2011. Addition¬ 
ally, ISchneider et al.l (1201511 find that the outer radial 
extent of the H 2 fluorescence luminosity doesn’t change 
significantly between epochs, which is a result consistent 
within the errors on our modeled v out estimates of the 
A A Tau epochs. 

The “extra absorber” obscures the inner disk H 2 flu¬ 
orescence in the 2013 HST-COS FUV spectrum, making 
AA Tau appear as a disk with a deficit of inner disk emis¬ 
sion - effectively, a pseudo-transition disk. Our modeling 
framework was capable of identifying the change in emis- 
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targets. Each radiation distribution was calculated from the minimum reduced-\' 2 model parameters that best reproduce the observed H 2 
emission lines. The top plot represents radiation distributions for all primordial disk targets, and the bottom plot shows the distributions 
for transition disk objects. The two disk evolution types appear to show an evolving H 2 FUV radiation field; primordial disks generally 
radiating more inward in the disk, with the bulk of the radiation occurring within r < 1 AU, and transition disk H 2 radiation starting at 
larger radii (r ~ 0.1 AU) extending to larger radii (r ~ 10 AU). 

sion line profiles between the 2011 and 2013 AA Tau 
observations and found that the bulk of the 2013 AA 
Tau [1,4] radiation in the disk originated at larger radii 
than the 2011 H 2 fluorescence. We expect our models 
are therefore capable of distinguishing between H 2 fluo¬ 
rescence evolution in differing disk types. 


5. DISCUSSION 

We created 2D radiative transfer models to simulate 
observed HST-COS and -STIS FUV H 2 emission lines 
to understand where the majority of the radiation arises 
in PPDs. We use the best-fit model results to define the 
inner and outer radii of warm H 2 (rj n , r out ) and examine 
if and how the molecular distributions change as PPDs 
evolve. We compare r m and r out to other dust and molec¬ 
ular tracers that help describe the evolutionary state of 
the PPDs. Table [I] provides a detailed list of inner disk 
observables for each target, including dust cavity radius 
(cavity ) and inner disk CO radius (ii n ,co)- We also look 
at where the theoretical snow lines in the disks exist and 
how these radii relate to the H 2 disk emission. 


5.1. Radiation Distribution of Modeled H 2 Fluorescent 

Emission 

Figure [9] presents the normalized radial distributions of 
warm H 2 transitioning from excited state [1,4] for all tar¬ 
gets. We modeled 6 primordial disks (AA Tau, BP Tau, 
DF Tau A, HN Tau A, RECX-11, and RECX-15) and 
8 transition disks (CS Glia, DM Tau, GM Aur, LkCa 
15, SU Aur, TW Hya, UX Tau A, and V4046 Sgr) to 
compare the radial distribution of warm H 2 in the disk 
atmospheres as the dust disk evolves. The H 2 radial dis¬ 
tributions of the different dust disk stages appear qual¬ 
itatively different. The primordial disk population (top 
plot in Figure O generally starts radiating significantly 
in the very inner disk (r < 0.05 AU), and the radiation 
only extends out to a few AU, consistent with the sim¬ 
ple _esthnates_of theaverage H 2 emitting radius presented 
bv lFrance et al.l (I2012bl) . The generalized transition disk 
radiation behavior (bottom plot) starts further out in 
the disk (r ~ 0.1 AU) and extend significantly further 
out into the disk (r ~ 10 AU). These different behaviors 
suggest structural changes in any of all of the following: 
the spatial distributions of warm H 2 in populations [u, J]; 
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Figure 10. A comparison of observed [1,4] progression line pro¬ 
files of targets with inclination angles between 30° and 40°. The 
two broadest line profiles, BP Tau and HN Tau A, are in the pri¬ 
mordial phase. The two narrowest line profiles, DM Tau and UX 
Tau A, are in the transition phase. 


the degree of Lyct penetration into the disk by clearing 
H 2 from the inner disk atmosphere; or the evolution of 
the disk surface temperature distribution. This evolving 
radiation structure is also observable in the line profiles 
of the [1,4] progression, as seen in Figure [TO] As the 
PPDs in our survey evolve from primordial to transition 
disks, the majority of the observed H 2 emission migrates 
to larger radii. 

We compare estimates of r in and r out to investigate 
the idea that the radial distributions of fluorescing H 2 
migrate outward in the disks as PPDs evolve. Figure [FT] 
presents a comparison of r in and r out , which shows the 
annulus of H 2 emission extending further out into the 
disk as the inner disk radius moves outward. A line can 
be fit to represent the relationship between the inner and 
outer radiating disk radii for our survey targets: 

login (r out (H 2 )) = 0.79 log 10 (r in (J? 2 )) + 1-39, (13) 

where both logi 0 (r in (H 2 )) and logi 0 (r o „ t (H 2 )) are 
in units of AU, and the coefficients [1.39 ± 0.22, 
0.79 ± 0.21] are computed from a y 2 minimization 
(Xmin = 0-896) of a linear function between logi 0 (rj n (H 2 )) 
and logio(?' ou t(H 2 )). The Spearman rank correlation co¬ 
efficient between r, n and r out indicates a statistically sig¬ 
nificant correlation between the variables (p = 0.70) with 
a small probability that the sample is randomized (n = 
5.5 x 10 -3 ), providing additional evidence that support 
the migration of the radial H 2 emission as PPD warm 
dust dissipates from the inner disk. 

5.2. Comparison to Dust Evolution 

We compare results from our modeled H 2 [1,4] progres¬ 
sion radial distributions with dust disk evolution diagnos¬ 
tics to gain insight into how the molecular inner disk envi¬ 
ronment of PPDs changes as dust grains clear. We iden¬ 
tify PPD evolution using observed color-color changes 
in the near- to rnid-IR SED slope of the disk, which pro- 
vi des an esti mate of the degree of warm dust clearing (see 
lEsnaillat et al.l[2014f) . We interpret the slope of each tar¬ 


Figure 11. We present the relation between the estimated rj n 
and Tout quantities, determined from the best-fit modeled radiation 
distributions for all targets. The inner radial boundary (r i n ) is 
defined as the inner radius of H 2 radiation in the disk that marks 
where at least 90% of the total radiation is accounted for in the 
outer disk. Likewise, the outer radial boundary (r 0 ut) is defined 
as the outermost radius of H 2 radiation that encompasses 90% of 
the total amount of radiation accounted for in the inner disk. The 
blue diamonds with error bars represent each modeled r i n and r ou t, 
and the black dashed line represents a linear fit to the data. The 
Spearman rank correlation coefficient ([p, n] = [0.70, 5.5 x 10 3 ]) 
between the two radial quantities suggest a strong increasing trend 
between them, indicating that the whole emitting region is moving 
outward. 



Figure 12. Comparison of r j n and r ou t with a n observable dust 
evolution diagnostic ni 3_31 ([Furlan et al.l[20090 . In the top plot: 
Each blue triangle with error bars represents each target point in 
our survey. The black dashed line represents the best-fit linear cor¬ 
relation between rj n and ni 3 _ 31 . In the bottom plot: Each green 
triangle with error bars represents each target point in our survey. 
The black dashed line represents the best-fit linear correlation be¬ 
tween Tout and ni 3 _ 3 i. In both plots, a clear increasing trend is 
seen in the radial H 2 emission boundaries as the warm dust disk 
content evolves. 


get SED with the observable 7113-31 ( Furlan et a l. 200 9:: 


log(\ 31 F\ 31 ) - log(\ 13 Fx 13 ) 
log( A 31 ) - log( A 13 ) 


(14) 


which is dominated by longer wavelength continuum 
emission from the optically-thick dust in the disk and is 
sensitive to the degree of dust settling towards the disk 
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Table 5 

Disk Parameters from Results & Literature 


Target 

ni3_3i M a 

^in,H 2 

^ out,H2 

r in,CO 

F cavity 

T(H 2 ) 

ref. b 


(IO” 8 Mq yr- 1 ) 

(AU) 

(AU) 

(AU) 

(AU) 

K 



AA Tau 

-0.51 

1.5 

0.08 ± 0.01 

BP Tau 

-0.58 

2.9 

0.04 ± 0.01 

CS Cha 

2.89 

5.3 

0.23 ± 0.05 

DF Tau A 

-1.09 

17.7 

0.04 ± 0.01 

DM Tau 

1.30 

0.29 

0.11 ± 0.01 

GM Am 

1.76 

0.96 

0.10 ± 0.01 

HN Tau A 

-0.44 

0.13 

0.04 ± 0.01 

LkCa 15 

0.62 

0.31 

0.20 ± 0.04 

RECX-11 

-0.80° 

0.03 

0.07 ± 0.02 

RECX-15 

-0.20 c 

0.10 

0.05 ± 0.01 

SU Aur 

0.74 

0.45 

0.35 ± 0.12 

TW Hya 

0.20 c 

0.02 

0.38 ± 0.14 

UX Tau A 

1.83 

1.00 

0.25 ± 0.06 

V4046 Sgr 

0.32 c 

1.30 

0.11 ± 0.01 


3.47 ± 0.54 

0.10 


4000l?*° 0 

2,11 

0.87 ± 0.10 

0.03 


2000±300 

1,2,11 

21.88 ± 4.68 


40 

2500±400 

5 

1.26 ± 0.89 

0.10 


1500li°°° 

11 

2.19 ± 1.48 


3 

2000±500 

3,4,10 

7.59 ± 2.75 

0.20 

20 

3000±450 

1,2,3,4,11 

3.80 ± 0.20 



2500±750 


6.03 ± 2.45 

0.10 

46 

1500t“g° 

1,3,6,10,11 

3.98 ± 2.00 



3000)( (250 


2.63 ± 1.62 


7.5 ± 1.5 

2500±350 

12 

12.02 ± 3.47 



1500±“®° 


3.98 ± 1.0 

0 l+°- 2 
u -0.04 

4 

2000±f°° 

3,4,10 

12.03 ± 3.46 

0.30 

25 

1500±“g° 

3,6,11 

3.31 ± 1.82 


14 

2000±500 

8,9 


a All M values taken from llnglebv et aLl (120131) . 

b (1) lAke son et aD (120051) ; (2) [Andrews & W illiams| J2007| ); (3) | A ndrews et al.| (1201U); (4) ICal vet et aH (120051) ; 
(5]_[Espaillat et al. (2007a|); (6) [E spaillat et al. (2007b); (7) France et al. (2012b); (8) Raps oiT et al.l d 201 51k (9) 
IRosenfeld et al.l (120131 ); (TO) ISalvk et al.l (120091); ( 1 1) ISalvk et al. d2011al ); (12 ) IWoitke et al.l (120111 ) 
c For ni 3 _ 3 i fim values not listed in lFur lan et al. (20091). we use Equation 1141 to estimate the observable from known or 
modeled dust SED. 


midplane (ID’Alessio et al.1 12006fl . For many targets i n 
this work, ni 3 _ 3 i were available in lFurlan et al.l (I2009D . 
but for targets not included in the iFurlan et al.l (1200911 
survey, we calculate 7713-31 with known or modeled disk 
SEDs (for ex ample, an intrica te m odel of V4046 Sgr SED 
was found bv lRosenfeld et af]l2013l ). We interpret the re¬ 
sults of ni 3 _ 3 i as follows: if 7113-31 < 0 , the inner dust 
disk is optically thick, essentially a primordial disk; if 
ni 3 — 3 i > 0 , the disk dust is optically thin, indicative of 
dust clearing or settling a nd evidence for PPD evolution 
into the transition state (lLadal 119871 : iStrom et al.l 119891 : 
lAndre fc Montmerld 11994b Table [T] provides a list of 
71 i3_3i values for all targets in this survey. 

A comparison of the [1,4] emission boundaries (r,; n , 
Tout) to rii 3 _ 3 i is made in Figure [12] The top figure 
shows the relationship between r,; n and 7713 - 31 , and the 
bottom figure shows r out versus 7713 - 31 . The triangles in 
both plots represent each target in our survey, and the 
black dashed line in each plot shows the linear correla¬ 
tion between rversus ni 3 _ 3 i and T out versus 7113 - 31 . 
It is apparent that the molecular inner and outer disk 
emission radii show a positive correlation with the dust 
disk evolution: the Spearman rank correlation coefficient 
for Ti n versus 7113-31 is 0.72 (n = 4.0 x 10~ 3 ), and p = 
0.69 for T out versus 7113-31 (71 = 6.9 x 10 -3 ). Both cor¬ 
relation coefficients suggest a strong increasing trend in 
the radial outward migration of the FUV H 2 radiation 
as the warm dust disk evolves in the disk samples. The 
linear correlation between rand 7113-31 is expressed as: 


We note that, for all transition disks in this study, r out 


10' 6 

IO ' 7 

| 10 ° 
•s 


10- 9 

IO 10 

0.01 0.10 1.00 
r, N (H 2 ) [AU] 

Figure 13. Comparison of the modeled inne r H 2 emission radiu s 
to the mass accretion rate of the target (from [Inglebv et al.ll2013l) . 
The purple x-points represent all targets with mass accretion rates 
> 10 -9 M® yr —1 , while the red diamonds represent the RECX 
targets (which have low mass accretion rates for primordial PPD 
targets). The black line is a negative correlation fit through all the 
purple points, suggesting that the mass accretion rate decreases as 
Tin increases. Since the accretion luminosity, more specifically the 
stellar Lya flux produced by the accretion, is directly related to 
the H 2 emission observed, it is important to note that rj n is not 
necessarily correlated to the flux produced by the mass accretion 
rate. Instead, r i n is sensitive to the observed emission line width, 
which is independent of the stellar incident flux. 


H ^ 
M 








~ 7 ' ft..' X 






l°gl( }( r in{H 2 )) =(0.19 ± 0.07) X 7113-31 

- (1.05 ±0.08), ( } 

and the linear correlation between r out and 7713-31 is ex¬ 
pressed as: 

^og 10 (r ou t(H 2 )) =(0.25 ±0.06) x ni 3 _ 3 i 

± (0.52 ±0.07). ( ' 


is found to be within the dust gap radius. One interpre¬ 
tation of this result, paired with the correlation between 
Tin and 7713 - 31 , is that the H 2 FUV radiation observed 
from the inner PPD atmosphere lags behind the dust 
disk evolution. 

This does not automatically mean that the molecular 
content of the disk is clearing, and we need further evi¬ 
dence of evolution with other inner disk molecular trac- 
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ers bef ore we can make this distinction. iFrance et al.1 
(l2012bf) outlined the conditions needed in the H 2 disk at¬ 
mosphere to produce Lya-pumped H 2 fluorescence. The 
opacity of absorbing H 2 in ground-state [u,J] must be 
large, with excitation temperatures T exc > 1500K, and 
the mass accretion rate (M) onto the proto-star must be 
large enough to produce enough Lya photons to stim¬ 
ulate the molecules. The mass accretion rate implies 
there is a reservoir of material in the inner regions of 
PPDs that feeds onto the proto-star, and a decrease in 
M over time (e.g., iMuzerolle et all [200 0) strongly sug¬ 
gests that the inner disk material is being depleted. Fig¬ 
ure M shows the relationship between M and iy n (H 2 ), 
with purple points representing r in (H 2 ) and M for all 
targets except the RECX targets, which are represented 
at re d diamonds. All mass accretion rates are taken from 
llnglebv et al.1 (1201311 . Figure [13] shows a negative cor¬ 
relation between M and i'i n (H 2 ), with Spearman rank 
correlation [p,n] = [-0.80, 1.9 x 10~ 3 ] (not including the 
RECX targets), suggesting that the H 2 atmosphere may 
be physically thinned or in different ground-state popula¬ 
tions not suitable for Lya-pumping in the very inner disk 
regions of evolved PPDs. The outlier points in Figure [T3l 
RECX-11 and RECX-15, appear to have abnormally low 
mass accretio n rates given the evolutionary stage of the 
disks (llnglebv et al.l 120111) . and more targets of varying 
evolution may be needed to understand if this result is 
universal among a large sampling of PPDs. It is im¬ 
portant to note that iis primarily dervied from the 
observed line widths of H 2 emission profiles, so determi¬ 
nation of r i n is largely independent of the incident FUV 
flux. 

The link between M and ri„(H 2 ) suggests that the in¬ 
ner disk is clearing of material as the mass accretion 
rate declines. One explanation for this correlation is 
that the warm H 2 atmosphere dissipates with the small 
dust grains. Dust grains present in the disk atmospheres 
of primordial disks may give warm H 2 a formation site 
to replenish molecul es lost to photo-dissociation and 
stellar accretion (se e IAugasonl [l970: Ha bart et al.11200'S 
iFleming et, alJl201fll h As the dust grains clear out and 
settle towards the disk midplane or evaporate from the 
inner disks of evolving PPDs, the warm H 2 atmosphere 
no longer has a formation site to maintain the molecular 
reservoir. Via accretion and photo-dissociative processes 
with FUV continuum photons between 912 - 1120 A, the 
leftover warm H 2 will continue to disperse, even as the 
accretion flux decreases. This leaves an optically thin 
(lV(if 2 ) < 10 18 cm -2 ) path for stellar Lya to reach the 
warm H 2 material at larger disk radii (r > 3 AU). 

The migration of r oltt (H 2 ) with increasing ni 3 _ 3 i also 
suggests that neutral hydrogen (HI) is being cleared from 
the inner disks of transitional PPDs. Photo-excitation 
via stellar Lya drives the H 2 fluorescence observed in the 
disk atmospheres, and as the emitting H 2 is observed fur¬ 
ther out in the disk, there must be new paths open for 
stellar UV radiation to reach the outer disk material. In 
primordial disks, HI re-processes and scatters incident 
stellar Lya down into the inner disk dFogel et al.ll201 ll ) 
while H 2 self-shields the radiation from penetrating to 
the outer disk, preventing the stellar Lya from reaching 
the outer disk effectively. If H 2 and HI column densities 


in the inner disk become optically thin in transitional 
disks, more stellar Lya can irradiate molecular material 
in the outer disk and may explain the observed correla¬ 
tion between r ollt (H 2 ) and ni 3 _ 3 i. This suggests that HI 
clearing from the inner disk may happen over a similar 
timescale as the charac teristic dust dissipation (iWvattl 
1 2 0081 iRibas et al] [2014 1 and mass accretion quenching 
( Fedele et al.ll201 3). This inner-to-outer disk dissipation 
is in agreement with the UV switch model, which de¬ 
scribes the dispersal of inner disk gas cut off from the gas 
reservoir of the outer disk, due to selec tive ph otoeva pora¬ 
tion of material out to r ~ 5 - 10 AU (jClarke et al.ll200 it 
lAlexander et al.ll2006f) . Observations of other outer-disk 
molecules photo-excited by Lya radiation provide addi¬ 
tional evidence for the loss of HI in the inner disks oftran- 
sitional objects. For example, iSchindhelm et al.l •: 2012 a; 
observe FUV-CO fluorescence, also powered by stellar 
Lya-pumping, at T exc ~ 500K, in transitional phase ob¬ 
jects with an average emission radius R co ~ 1 - 5 AU. 
This indicates that less HI and H 2 column is present in 
the inner disk to shield the stellar Lya flux from reaching 
the cooler CO material at intermediate radii in transition 
systems. 

Figure [14] shows a ID radial comparison of dust and 
molecular tracers determined in our targets. We present 
the locations of the outer radiation boundary for H 2 
FUV emission, as determined from our models (r outj [-i, 4 i; 
green triangles), and the observed dust cavity walls of the 
transitional disk populations (jcavity ; blue squares). For 
all transitional disks, we find r oztt; [ 1 , 4 ] inward of r cav ity, 
meaning that the H 2 population observed in all transi¬ 
tion PPDs radiates where the dust is optically thin, sug¬ 
gesting that the H 2 populations remain optically thick 
even after the dust gr ains have dissipated. Studies like 
Ivan der Marel et aD (12015 ) also find a substantial deple¬ 
tion of the dust-to-gas ratio inside the dust cavities of 
well-studied transition disks, which is consistent with our 
findings. 

5.3. Near-IR CO Emission and Comparison to Snow 
Line Radii 

Figure [14] includes radial estimates of the inner 
radiation boundary for H 2 FUV emission (r irl) [ 1 , 4 ] 5 
blue x’s) and the inner radius of near-IR CO emission, 
determ ined from LTE models presented by iSalvk et al.1 
(201Ta) (r in,co', black diamonds). The inner disk 
emission radii of FUV H 2 and near-IR CO appear to 
be^ roughly co -spat ial, which is a result also found by 
IFrance et al.l (|2012bD when comparing the observed 
FWHMs of FUV H 2 fluorescence e mission an d n ear-I R 
CO emission. An extensive study bv lBrown et al. (20131 ! 
concluded there is a correlation between the near-IR 
CO P(8) equivalent width and dust disk dispersal in 
transitional disks, suggestive of outer radial origins of 
the CO emission as PPD dust evolves. We have shown 
that i'i n ,[i, 4 ] increases with ni 3 _ 3 i and decreases with 
M, providing further evidence that the inner gas disk 
environment becomes optically thin as disks evolve 
towards the transition stage. 

We note the disk locations of possible theoretical 
snow lines in PPDs and these radii coincide with the 
H 2 fluorescence in Figure [T4j As the disk evolves, it 
cools over time, so the snow line is expected to migrate 






















































The Evolution of Gas in Transition Disks 


15 


TW Hya 
LkCal5 
V4046 SGR 
UX Tau A 
SU Aur 
GM Aur 
DM Tau 
CS Cha 
RECX-15 
RECX-11 
HN Tau 
DF Tau 
BP Tau 
A A Tau 


0.01 0.10 1.00 10.00 

Gisk (ATJ) 


♦ Geo ♦ ■ X ■ 



■ bavity ♦ 


* ■ 

X r in ,[i,4] $$ 

■ 

^ Git, [1,4] 


a m 

+— H- 

V//A 

A 

* ♦ 

a m 

« 

■ A ■ 



a m 

X 


m - 

H 



X 


A 

X ♦ — 

-*—■ 

— 

♦X * 

\///\ 

— 





Figure 14. A radial comp arison of the inner and outer extent of FUV H 2 emission (this work), the innermost radius of near-IR CO 
emission (|Salvk et al.N2011afl . and dust cavity locations in transition disk targets (see Table [T| for references). The light blue shaded area 
from 1 AU < r dis k < 3 AU represen ts the theoretical water-ice snow line for the presence of water-ice at the midplane of primordial and 
transitional PPDs (jBaillie et al.ll2015l) . 


inward in the disk as the protostellar syst em age s (?). 
Several independent studies (e.g. iMeiierink et al.l 120091 : 
iMandell et al.l I2012H conclude that the location of the 
water-ice snow line in PPDs are expected to be found 
within r ~ 1 - 3 AU for all PPD states. iBaillie et al.l 
(I20l5h shows that the evolution of the water-ice snow 
line at all stages of PPD evolution (from ages 10 6 - 
f0 7 yr) only varies by ~ 0.5 AU. Observations of H 2 O 
and OH (which is thought to be a bi-product of H 2 0 
photo-dissociation) in the near- and mid-IR are also 
consi s tent with these condensation radii (IMalfait et al. 
1998 ; iCarr et all 12004 IMandell et al.l 120081 : iSalvk et al. 
20081 1. Figure [lj] includes a shaded blue region that 
represents the assumed generalized H 2 0 snow line radii 
in PPDs, located between r disk = 1-3 AU. With 
the exception of BP Tau and DF Tau A, all targets 
have outer H 2 emission radii that extend to within or 
outward of the water-ice snow line. 

6 . CONCLUSIONS AND FUTURE WORK 

We have created 2D radiative transfer models of FUV 
H 2 fluorescence emission in PPDs and compared them 
with observations made with HST -COS and STIS. We 


analyze the radial distribution of H 2 emission produced 
by parametrized models, which are determined using a 
reduced-x 2 statistic, to understand how the emitting H 2 
regions changes as PPD dust disks evolve. We summarize 
our findings and interpret the evolutionary behavior of 
the molecular disk atmosphere as the inner dust disks of 
PPDs disperse: 

1. The modeled H 2 radial distributions differ between 
primordial and transitional disks. Primordial disks 
have the majority of the total H 2 flux arising from 
the innermost disk radii and less produced outside 
~ 1 AU. For transitional disks, the total H 2 flux 
migrates to larger disk radii, producing less flux in 
the innermost disk and more out to r ~ 10 AU. 

2. We see a positive correlation between the result¬ 
ing inner and outer emission radii of FUV H 2 (rj n 
and Tout), which supports the result described in 
conclusion 1. This can mean: a) that the physical 
structure (i.e., temperature) of the warm molecular 
disk atmosphere changes as PPDs evolve, b) the 
warm, ground-state H 2 populations [u,J] change, 
resulting in evolving regions of the disks where the 
warm H 2 atmosphere will reprocess the stellar Lyct 
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radiation field, or c) H 2 is being destroyed in the 
inner disk and not re-formed, owing to the lack of 
dust grains; the latter point allows stellar Lya to 
penetrate to larger r out . 

3. We observe positive correlations between r ln , r out , 
and ni 3 _ 3 i, suggesting that rj„ corresponds with 
the loss of warm, small dust grains in the inner¬ 
most disk. We find a negative correlation between 
ri n and M, providing evidence that the warm H 2 
inner disk atmosphere may be physically thinned 
or cleared as the PPDs evolve, possibly by the loss 
of a molecular formation site as the dust grains 
dissipate from the atmosphere. Using the observed 
dust cavity radii of the transitional disk targets, we 
compare r ollt to r cav it y and find that, for all tran¬ 
sition disk targets, r oui is found inward of r cavity- 
This indicates that the warm H 2 disk (for r > r i n ) 
remains optically-thick where the warm dust grains 
are optically-thin in the disks. This suggests that 
the physical mechanism that clears or settles the 
inner disk dust either does not have the same ef¬ 
fect on the molecular disk atmosphere, or there is a 
time lag for the gas disk to respond to the changes 
observed in the dust distribution. 

4. We examine where the emitting H 2 originates in 
the disk relative to warm CO and the theoretical 
location of water-ice snow lines. Inner disk CO is 
roughly co-spatial with iy n for all targets, which 
could point to the dispersal of the warm molecular 
disk atmospheres of evolving disk systems. With 
the exception of a few primordial disk targets, all 
targets have emitting H 2 regions that encapsulate 
the theoretical water-ice snow line. If disk clear¬ 
ing mechanisms, such as disk photoevaporation via 
EUV/X-ray photons, are primarily responsible for 
the final dispersal of the gas disk at the end of the 
PPD lifetime, it is important to examine late-type 
PPDs to monitor molecular disk clearing as transi¬ 
tional disks evolve to debris disks. 
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